extern doublereal scnrm2_(integer *n, complex *x, integer *incx);
